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We consider a quantum system consisting of a regular chain of elementary subsystems with nearest 
neighbor interactions and assume that the total system is in a canonical state with temperature T. 
We analyze under what condition the state factors into a product of canonical density matrices 
with respect to groups of n subsystems each, and when these groups have the same temperature T. 
While in classical mechanics the validity of this procedure only depends on the size of the groups 
n, in quantum mechanics the minimum group size n m i n also depends on the temperature T ! As 
examples, we apply our analysis to a harmonic chain and different types of Ising spin chains. We 
discuss various features that show up due to the characteristics of the models considered. For the 
harmonic chain, which successfully describes thermal properties of insulating solids, our approach 
gives a first quantitative estimate of the minimal length scale on which temperature can exist: This 
length scale is found to be constant for temperatures above the Debye temperature and proportional 
to T~ 3 below. 
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I. INTRODUCTION 

Thermodynamics is among the most successfully 
and extensively applied theoretical concepts in physics. 
Notwithstanding, the various limits of its applicability 
are not fully understood Q, Q ■ 

Of particular interest is its microscopic limit. Down to 
which length scales can its standard concepts meaning- 
fully be defined and employed? 

Besides its general importance, this question has 
become increasingly relevant recently since amazing 
progress in the synthesis and processing of materials with 
structures on nanometer length scales has created a de- 
mand for better understanding of thermal properties of 
nanoscale devices, individual nanostructures and nanos- 
tructured materials 0. Q- 0- Experimental techniques 
have improved to such an extent that the measurement of 
thermodynamic quantities like temperature with a spa- 
tial resolution on the nanometer scale seems within reach 




To provide a basis for the interpretation of present day 
and future experiments in nanoscale physics and tech- 
nology and to obtain a better understanding of the lim- 
its of thermodynamics, it is thus indispensable to clarify 
the applicability of thermodynamical concepts on small 
length scales starting from the most fundamental theory 
at hand, i. e. quantum mechanics. In this context, one 
question appears to be particularly important and in- 
teresting: Can temperature be meaningfully defined on 
nanometer length scales? 

The existence of thermodynamical quantities, i. e. the 
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existence of the thermodynamic limit strongly depends 
on the correlations between the considered parts of a sys- 
tem. 

With increasing size, the volume of a region in space 
grows faster than its surface. Thus effective interactions 
between two regions, provided they are short ranged, be- 
come less relevant as the sizes of the regions increase. 
This scaling behavior is used to show that correlations 
between a region and its environment become negligible 
in the limit of infinite region size and that therefore the 
thermodynamic limit exists 0, ITTL IT^| . 

To explore the minimal region size needed for the appli- 
cation of thermodynamical concepts, situations far away 
from the thermodynamic limit should be analyzed. On 
the other hand, effective correlations between the consid- 
ered parts need to be small enough 0,0]. 

The scaling of interactions between parts of a system 
compared to the energy contained in the parts themselves 
thus sets a minimal length scale on which correlations 
are still small enough to permit the definition of local 
temperatures. It is the aim of this paper to study this 
connection quantitatively 

Some attempts to generalize thermodynamics such 
that it applies to small systems have been made 0, Q, 
Il5j . These approaches consider ensembles of indepen- 
dent, i.e. nonintcracting, small systems. By introducing 
an additional thermodynamical potential they take into 
account the surface effects of the small systems. However, 
since the interactions between the small systems are ne- 
glected, these concepts cannot capture the physics of the 
correlations. This shortcoming is also obvious from the 
results: The correction terms they predict do not depend 
on temperature, whereas it is well known, that correla- 
tions become more important the lower the temperature. 

Recently the impact of quantum correlations, i. e. en- 
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tanglcmcnt on macroscopic properties of solids and phase 
transitions has drawn considerable attention |2(J, |21j . 
Since our analysis of criteria for local temperatures is 
based on a study of correlations, our theoretical approach 
is a promising tool to provide further insight into the role 
of correlations in solid state physics. 

We adopt here the convention that a local tempera- 
ture exists if the considered part of the system is in a 
canonical state, where the distribution is an exponen- 
tially decaying function of energy characterized by one 
single parameter. This implies that there is a one-to- 
one mapping between temperature and the expectation 
values of observables, by which temperature is usually 
measured. Temperature measurements based on differ- 
ent observables will thus yield the same result, contrary 
to distributions with several parameters. In large sys- 
tems composed of very many subsystems, the density of 
states is a strongly growing function of energy |22j | . If the 
distribution were not exponentially decaying, the prod- 
uct of the density of states times the distribution would 
not have a pronounced peak and thus physical quantities 
like energy would not have "sharp" values. 

There have been atempts to describe systems which 
are not in an equilibrium state but in some sense close 
to it with a generalized form of thermodynamics, that 
has additional system parameters. Such a situation ap- 
pears for example in glasses [T^ | . Our approach analyzes 
whether thermodynamics in its standard form can apply 
locally. A study whether a generalized form of thermody- 
namics might apply even more locally should be a subject 
of future research. 

A typical setup where the minimal length scale we cal- 
culate becomes relevant could be the measurement of a 
temperature profile with very high resolution etc. One 
is thus interested in scenarios where the entire sample is 
expected to be in a stationary state. In most cases this 
state is close to a thermal equilibrium state [23J. 

Based on the above arguments and noting that a quan- 
tum description becomes imperative at nanoscopic scales, 
the following approach appears to be reasonable: Con- 
sider a large homogeneous quantum system, brought into 
a thermal state via interaction with its environment, di- 
vide this system into subgroups and analyze for what 
subgroup-size the concept of temperature is still applica- 
ble. 

Harmonic lattice models are a standard tool for the 
description of thermal properties of solids. We therefore 
apply our theory to a harmonic chain model to get esti- 
mates that are expected to be relevant for real materials 
and might be tested by experiments. 

Recently, spin chains have been subject of extensive 
studies in condensed matter physics and quantum in- 
formation theory. Thus correlations and possible local 
temperatures in spin chains are of interest, both from a 
theoretical and experimental point of view [2~iL [25| . We 
study spin chains with respect to our present purpose and 
compare their characteristics with the harmonic chain. 

This paper is organized as follows: In section [HJ we 



present the general theoretical approach which derives 
two conditions on the effective group interactions and 
the global temperature. In the following two sections we 
apply the general consideration to two concrete models 
and derive estimates for the minimal subgroup size. Sec- 
tion II I II deals with a harmonic chain, a model with an 
infinite energy spectrum. In contrast, a spin chain has a 
bounded energy spectrum. Section llVl thcrcfore discusses 
an Ising spin chain in a transverse field. In the conclu- 
sions section El we compare the results for the different 
models considered and indicate further interesting topics. 

II. GENERAL THEORY 

We consider a homogeneous (i.e. translation invariant) 
chain of elementary quantum subsystems with nearest 
neighbor interactions. The Hamiltonian of our system is 
thus of the form , 

H = J2 H i+ki+i (!) 

i 

where the index i labels the elementary subsystems. Hi 
is the Hamiltonian of subsystem i and the interac- 

tion between subsystem i and i + 1. We assume periodic 
boundary conditions. 

We now form N@ groups of n subsystems each (index 
i -> (/i - l)n + j; n = 1, . . . , N G ; j = 1, . . . , n) and split 
this Hamiltonian into two parts, 

H = H a + I, (2) 

where Hq is the sum of the Hamiltonians of the isolated 
groups, 

Hq = i^Hfi — I(m,(m+l) with 

n 

W M = ^ -Hn(jU-l)+3 + ^w(M-l)+j,w(M-l)+j+l (3) 

and / contains the interaction terms of each group with 
its neighbor group, 

N G 

1 = X! ( 4 ) 

We label the eigenstates of the total Hamiltonian H and 
their energies with the Greek indices (if, ip) and eigen- 
states and energies of the group Hamiltonian Ho with 
Latin indices (a, b), 

H\<p)=E lp \<p) and H \a)=E a \a). (5) 

Here, the states |a) are products of group eigenstates 

N G 

\a) = II K>> (6) 

where (H^ - i^/m+i) \a-n) = is the energy 

of one subgroup only and E a = ^n- 
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A. Thermal State in the Product Basis 

We assume that the total system is in a thermal state 
with the density matrix 



(7) 



in the eigenbasis of H . Here, Z is the partition sum and 
(3 = (fcfiT) -1 the inverse temperature with Boltzmann's 
constant ks and temperature T. Transforming the den- 
sity matrix 10 into the eigenbasis of Hq we obtain 



(a\p\a) 



W a (E) 



-pE 



Z 



dE 



(8) 



for the diagonal elements in the new basis. Here, the 
sum over all states \tp) has been replaced by an integral 
over the energy. E is the energy of the ground state 
and Ei the upper limit of the spectrum. For systems 
with an energy spectrum that does not have an upper 
bound, the limit E\ — > oo should be taken. The density 
of conditional probabilities w a {E) is given by 



W a {E) 



1 

~ae 



E N^l 2 ( 9 ) 

{\ V ):E<E V <E+AE} 



where AE is small and the sum runs over all states \ip) 
with eigenvalues E v in the interval [E,E + AE). To 
compute the integral of equation © we need to know 
the distribution of the conditional probabilities w a {E). 

The state \a) is not an eigenstate of the total Hamilto- 
nian H. Thus, if H would be measured in the state |a), 
eigenvalues of H would be obtained with certain proba- 
bilities: w a (E) is the density of this probability distribu- 
tion. Since the hamiltonian H is the sum of hamiltonians 
of the groups, the situation has some analogies to a sum 
of random variables. This indicates that there might ex- 
ist a central limit theorem for the present quantum sys- 
tem, provided the number of groups becomes very large 
|28j . Since the state \a) is not translation invariant and 
since H also contains the group interactions, the central 
limit theorem has to be of a Lyapunov (or Lindeberg) 
type for mixing sequences [27||. One can indeed show 
that such a quantum central limit theorem exists for the 
present model [2^, |^| and that w a (E) thus converges 
to a Gaussian normal distribution in the limit of infinite 
number of groups Nq, 



r 1 \ (E-E a -e a ) 
lim w a {E) = —== cxp — 



where the quantities e a and A a are defined by 

e a = (a\H\a) - E a 
A 2 = (a\H 2 \a) -(a\H\a) 2 . 



(10) 

(11) 
(12) 



e a is the difference between the energy expectation value 
of the distribution w a (E) and the energy E a , while is 



the variance of the energy E for the distribution w a (E) . 
Note that e a has a classical counterpart while A^ is 
purely quantum mechanical. It appears because the com- 
mutator [H 1 Hq] is nonzero, and the distribution w a (E) 
therefore has nonzero width. The two quantities e a and 
A^ can also be expressed in terms of the interaction only 
(see eq. ©), 



e a = (a\I\a) 
Al = 



a\I 2 \a) - (a\I\a)\ 



(13) 
(14) 



meaning that e a is the expectation value and A 2 the 
squared width of the interactions in the state |a). 

The rigorous proof of equation ljTU|l is given in [2^ and 
based on the following two assumptions: The energy of 
each group as defined in equation © is bounded, i. 
e. 

(X\HM<C (15) 
for all normalized states |x) and some constant C, and 
(a\H 2 \a) - (a\H\a) 2 > N G C (16) 

for some constant C > 0. 

In scenarios where the energy spectrum of each ele- 
mentary subsystem has an upper limit, such as spins, 
condition (|15|) is met a priori. For subsystems with an 
infinite energy spectrum, such as harmonic oscillators, 
we restrict our analysis to states where the energy of ev- 
ery group, including the interactions with its neighbors, 
is bounded. Thus, our considerations do not apply to 
product states \a), for which all the energy was located 
in only one group or only a small number of groups. Since 
Nq 3> 1, the number of such states is vanishingly small 
compared to the number of all product states. 

If conditions (|15fl and (|16|1 are met, equation JHJ can 
be computed for Nq ^> 1 |30j : 



1 / (i 2 A 2 
( a \p\ a ) = exp -(3y a H —2- 



crfc 



2Z 

E -y a + pA 2 a 
V2A a 



erfc 



Eh, -y a + pA 2 a 
V2A a 



(17) 

where y a = E a + e a and erfc(a;) is the conjugate Gaussian 
error function, 



erfc(x) = —= 



ds. 



(18) 



The second error function in (|17fl only appears if the 
energy is bounded and the integration extends from the 
energy of the ground state Eq to the upper limit of the 
spectrum E\. 

Note that y a is a sum of Nq terms and that A a ful- 
fills equation (|16|l . The arguments of the conjugate error 
functions thus grow proportional to s/Nc or stronger. 
If these arguments divided by \f Ng arc finite (different 
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from zero), the asymptotic expansion of the error func- 
tion may thus be used for N G S> 1: 



erfc(x) 



cxp ( — x 



'■K X 



for ; 



x — ► 00 



exp i—x 2 ) 
2 + — ^= '- for x -> -00 

s/TTX 



Inserting this approximation into equation (|17|l and using 
Eo < y a < Ei shows that the second conjugate error 
function, which contains the upper limit of the energy 
spectrum, can always be neglected compared to the first, 
which contains the ground state energy. 

The same type of arguments show that the normaliza- 
tions of the Gaussian in equation IjlUfl is correct although 
the energy range does not extend over the entire real axis 
(—00, 00). 

Applying the asymptotic expansion (|19fl , equation l|17|) 
can be taken to read 



(a\p\a) 



Z 



cxp 



P E a 



(20) 



for (Eo -E a -e a + 0A 2 a ) / (V27Vg A„) < and 



exp 



(a\p\a) 



~0E O 



(E a + e a - E ) 2 
2A 2 



2ttZ 



E a -E a -e a + (3Aj 

Aa 



(21) 



for (Eo -E a -e a + (3A 2 a ) / (^2N^ A a ) > 0. 

The off diagonal elements (a\p\b) vanish for 
I E a — Eb I > A a + Aft because the overlap of the two distri- 
butions of conditional probabilities becomes negligible. 
For \E a — Eb\ < A a + At, the transformation involves 
an integral over frequencies and thus these terms are 
significantly smaller than the entries on the diagonal. 



B. Conditions for Local Thermal States 

We now test under what conditions the density matrix 
p may be approximated by a product of canonical den- 
sity matrices with temperature (3i oc for each subgroup 
p = 1,2,..., Nq. Since the trace of a matrix is invariant 
under basis transformations, it is sufficient to verify the 
correct energy dependence of the product density matrix. 
If we assume periodic boundary conditions, all reduced 
density matrices are equal and their product is of the 
form (o||3|a) oc exp(— /3i oc E a ). We thus have to verify 
whether the logarithm of rhs of equations (|20|l and l|21|) 
is a linear function of the energy E a , 



In ((a\p\a)) « -/3; oc E a + c, 



(22) 



where (3i oc and c are constants. 

Note that equation IL'l'l) does not imply that the occu- 
pation probability of an eigenstate \tp) with energy E v 
and a product state with the same energy E a w E v 



are equal. Since /3/ oc and f3 enter into the exponents 
of the respective canonical distributions, the difference 
between both has significant consequences for the occu- 
pation probabilities; even if /?/ oc and (3 are equal with 
very high accuracy, but not exactly the same, occupation 
probabilities may differ by several orders of magnitude, 
provided that the energy range is large enough. 

We exclude negative temperatures (/3 > 0). Equation 
(|22|l can only be true for 



E a + £a — Eq 



(23) 



as can be seen from equations Q2(Jfl and (|21|) . In this 
case, (a\p\a) is given by (|20|l and to satisfy l|22|) . e a and 
A^ furthermore have to be of the following form: 



?A 2 
2 a 



c\E a + c 2 



(24) 



where ci and C2 are constants. Note that e a and A 2 need 
not be functions of E a and therefore in general cannot be 
expanded in a Taylor series. 

To ensure that the density matrix of each subgroup p 
is approximately canonical, one needs to satisfy (|24fl for 
each subgroup p separately; 



(3 x: 



Cl En 



01 



where e M = (al^n^n+ik) with e a = Y,^=i £ ^^ 



(25) 



,a\K\a) (a\H^\a) 2 and A^ 



A 

EZZLMnu-in u + HvH v -i\a)-2(a\H u -i\a){a\H u \a 
and 

Temperature becomes intensive, if the constant C\ van- 
ishes, 



Icil « 1 



Ploc = P. 



(26) 



If this was not the case, temperature would not be inten- 
sive, although it might exist locally. 

It is sufficient to satisfy conditions (|23|l and (|25|l for 
an adequate energy range E min < E^< E max only. For 
large systems with a modular structure, i.e. a system 
composed of a large number of subsystems, the density 
of states is typically a rapidly growing function of energy 
[22Ll33T |. If the total system is in a thermal state, occupa- 
tion probabilities decay exponentially with energy The 
product of these two functions is thus sharply peaked at 
the expectation value of the energy E of the total system 
E + E =Tr(Hp), with E being the ground state energy 
(see figure P>. The energy range thus needs to be cen- 
tered around this peak and large enough. On the other 
hand it must not be larger than the range of values E^ 
can take on. Therefore a pertinent and "safe" choice for 
Emin and E max is 



E m in = max 
E m ax = min 



1 1 E 

Plmin ' a Nq 

E 



t^hnax 1 



<V<; 



Eq \ 

' N G J 

Eg \ 

N G J 



(27) 
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FIG. 1: The product of the density of states r](E) times the 
occupation probabilities (ip\p\ip) forms a strongly pronounced 
peak at E = E. 



where a 3> 1 and E will in general depend on the global 
temperature. In equation (|77|) . [E^] min and de- 
note the minimal and maximal values E^ can take on. 

Figure |21 shows the logarithm of equation (|17fl and the 
logarithm of a canonical distribution with the same (3 
for a harmonic chain as an example. The actual density 
matrix is more mixed than the canonical one. In the 
interval between the two vertical lines, both criteria 123|) 
and (|25|l are satisfied. For E < Ei ow 1(23(1 is violated and 
((25(1 for E > Ehigh- To allow for a description by means 
of canonical density matrices, the group size needs to be 
chosen such that E low < E min and E high > E max . 



\n((a\p\a)) 




E 



FIG. 2: In ((a|/3|a)) for p as in equation 1171 (solid line) and a 
canonical density matrix p (dashed line) for a harmonic chain. 

For a model obeying equations (|15|l and (|16|l . the two 
conditions 12-il) and (|25|) . which constitute the general 
result of this article, must both be satisfied. These fun- 
damental criteria will now be applied to some concrete 
examples. 



III. HARMONIC CHAIN 

As a representative for the class of systems with an 
infinite energy spectrum, we consider a harmonic chain 
of Ng ■ n particles of mass m and spring constant \pm loq. 



In this case, the Hamiltonian reads 

Hi = y^ + y^ (28) 

I iyi+1 = -mw 2 4 9Hi, (29) 

where pi is the momentum of the particle at site i and 
qi the displacement from its equilibrium position i ■ clq 
with ao being the distance between neighboring parti- 
cles at equilibrium. We divide the chain into Nq groups 
of n particles each and thus get a partition of the type 
considered above. 

The Hamiltonian of one group is diagonalizcd by a 
Fourier transform and the definition of creation and an- 
nihilation operators al and for the Fourier modes (see 
appendix . 

Ng , ,x 

£« = EE^ K(m) + 2 ). ( 3 °) 

|U=1 k ^ ' 

where k = Trl/(a (n +1)) (I = 1,2,..., n) and the fre- 
quencies Uk are given by io\ = 4wq sin 2 (fca/2). is 
the occupation number of mode k of group fi in the state 
| a). We chose units, where h = l. 

We first verify that the harmonic chain model fulfills 
the conditions for the applicability of the quantum cen- 
tral limit theorem I jlOj l. To see that it satisfies the con- 
dition l|16fl one needs to express the group interaction 
Ifm,iMn+i in terms of a k and afc, which yields A M = for 
all fj, and therefore 

N G 

A^ = £ A '' (31) 
where A M , the width of one group interaction, reads 




A?, has a minimum value since all rih > and all m v > 0. 
In equation yzj), k labels the modes of group p: with 
occupation numbers and p the modes of group \i + 1 
with occupation numbers m p . The width A 2 thus fulfills 
condition 1)16(1. 

Since the spectrum of every single oscillator is infinite, 
condition (|15fl can only be satisfied for states, for which 
the energy of the system is distributed among a substan- 
tial fraction of the groups, as discussed in section ITT1 

We now turn to analyze the two criteria Il2.'i[l and 1(25(1 . 
The expectation values of the group interactions vanish, 
Efi = 0, while the widths A 2 depend on the occupation 
numbers nj. and therefore on the energies E^. We thus 
have to consider both conditions, 1(23(1 and l(25|) . To an- 
alyze these, we make use of the continuum or Dcbyc 
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approximation |34| . requiring n > 1. oo < I, where 
I = noo, and the length of the chain to be finite. In 
this case we have u>k = V k with the constant velocity of 
sound v = LUoao and cos(fcao/2) as 1. The width of the 
group interaction thus translates into 



(33) 



where n + 1 ~ n has been used. The relevant energy 
scale is introduced by the thermal expectation value of 
the entire chain 



E = N G nk B G 



e/T 



■ dx, 



and the ground state energy is given by 
2 f e/r 



( | 



(i 



N G nk B <d 



(34) 



(35) 



We first consider the criterion 1)23)1 . 

For a given E a = Ylu E ^' the squared width A 2 is 
largest if all are equal, E^ = E V/i. Thus is 
hardest to satisfy for that case, where it reduces to 



E 



E^ 
N G 



> o. 



(36) 



Equation 1)36)) sets a lower bound on n. For temperatures 
where E < Eo, this bound is strongest for low energies 
E, while at E > Eo it is strongest for high energies E. 
Since condition 1)25)1 is a stronger criterion than condition 
(|23|l for E > Eo, we only consider 1)36)1 at temperatures 
where E < Eo- In this range, 1)36)1 is hardest to satisfy 
for low energies, i.e. at E = (E/aN G ) + (Eo/N G ), where 
it reduces to 



n > 



9 a 
T 4g 



4e 

Q 



(37) 



with e = E/(nN G k B <d)- 

To test condition 1)2 oil we take the derivative with re- 
spect to Efj, on both sides, 

where we have separated the energy dependent and the 
constant part in the lhs. (|38|l is satisfied if the energy 
dependent part is much smaller than one, 







^ [ E^x + E, +1 -2j^)<6<1 



(39) 



This condition is hardest to satisfy for high energies. Tak- 
ing E^-i and £^+i equal to the upper bound in equation 
lf27)l. it yields 



2a 9_ 



(40) 



where the "accuracy" parameter S <C 1 quantifies the 
value of the energy dependent part in 138(1 . 

Since the constant part in the lhs of (|38fl satisfies 



VJfi _Vs) 



n 2 N, 



a 



V2 



a 



(41) 



temperature is intensive. 

Inserting equation l|34|) into equation (|37|l and (|40|) one 

can now calculate the minimal n for given 5, a, and T. 



Figure |3] shows 



for a — 10 and S = 0.01 given by 



criterion l|37|) and l|40|l as a function of T/0. Hence, 
local temperature exists, i. e. local states are canonical 
for all group sizes larger then the maximum of the two 
^min-curves plotted in figure [3] 





10 2 



10 1 10 2 10 3 



T/e 



FIG. 3: Log-log-plot of n m i n from eq. 137H (dashed line) and 
Timin from eq. Il 101 (solid line) for a — 10 and 8 = 0.01 as a 
function of T/Q for a harmonic chain. S and a are defined 
in equations 1401 and 1271 . respectively. Local temperature 
exists in the shaded area. 



For high (low) temperatures n Tl 
mated by 



can thus be esti- 



2 a/ '6 
(3a/27r 2 ) 



for T > 9 
for T < 9 



(42) 



Equation (|42|) also shows the dependence of the results 
on the "accuracy parameters" a and S. In the whole tem- 
perature range, n m i n cx a, in other words, the larger one 
chooses the energy range where Q23JI and 1)25)1 should be 
fulfilled, the larger has to be the number of particles per 
group. Furthermore, for high temperatures, n m i n oc 
which simply states that one needs more particles per 
group to obtain a canonical state with better accuracy. 

Since the resulting minimal group sizes n m i n are larger 
than 10 3 for all temperatures, the application of the De- 
bye approximation is well justified. 
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IV. ISING SPIN CHAIN IN A TRANSVERSE 
FIELD 

In this section we consider an Ising spin chain in a 
transverse field. For this model the Hamiltonian reads 



Hi 



H,i+1 



-Bat 

J X 



y a". 



'i+i 



(43) 



where erf, a\ and of arc the Pauli matrices. B is the mag- 
netic field and J x and J v are two coupling parameters. 
We will always assume B > 0. 

The entire chain with periodic boundary conditions 
may be diagonalized via successive Jordan- Wigncr. 
Fourier and Bogoliubov transformations (see appendix 
IBl . The relevant energy scale is introduced via the ther- 
mal expectation value (without the ground state energy) 



-g _ nN G 
2ir 



dk 



where u> k is given in equation {B 
energy Eq is given by 

nN G 



The ground state 



2tt 



dk 



(45) 



Since N G 3> 1, the sums over all modes have been re- 
placed by integrals. 

If one partitions the chain into N G groups of n sub- 
systems each, the groups may also be diagonalized via 
a Jordan- Wigner and a Fourier transformation (see ap- 
pendix [Bj|. Using the abbreviations 



K 



Jx ~\~ J y 



and L 



Jx Jy 

2B 



(46) 



IB 

the energy E a reads 

No / 1\ 

E ° = 2B E E t 1 - K cos ( fe )] ( - o ) . ( 4? ) 

where k = nl/(n + 1) (Z = 1, 2, . . . , n) and n^(/Lt) is the 
fermionic occupation number of mode k of group /i in the 
state | a). It can take on the values and 1. 

For the Ising model at hand one has, as for the har- 
monic chain, e a = for all states |o), while the squared 
variance A 2 reads 



N G 

a 2 = Va 2 , 

Jt=l 



with 



-2B 2 (if 2 - L 2 



5 2 
2 



2 



L 2 
2 



(48) 



(49) 



T £sin 2 (fc) 



T ^sin 2 (p) (n£(/i+l) 



where the n£ (/i) are the same fermionic occupation num- 
bers as in equation (|4TJl . 

The conditions for the central limit theorem are met 
for the Ising chain apart from two exceptions: Condition 
(|15|l is always fulfilled as the Hamiltonian of a single spin 
has finite dimension. As follows from equation 1)49)1 . con- 
dition 1)16)1 is satisfied except for one single state in the 



case where J x = J v (L = 0) and J x 



-Jy (K = 0) 



respectively. These two states have A 2 = and thus 
A 2 < NgC". The state for L = is the one where 
all occupation numbers n%(fJ>) vanish and the state for 
K = is the state with alternating occupation numbers 
n%(ji) = 0, n%{ii + 1) = 1, n%(n + 2) = 0, . . . (for all 
k each). As there is, for given parameters, at most one 
state that docs not fulfill 1)16)1 , the fraction of states where 
our theory does not apply is negligible for N G ^> 1. 

We now turn to analyze conditions 1)23)1 and (|25(l . Since 
the spectrum of the Ising chain is limited, there is no ap- 
proximation analog to the Debye approximation for the 
harmonic chain and A 2 cannot be expressed in terms of 
£7 M _i and E^. We therefore approximate (|23fl and 1)25)1 
with simpler expressions. The results are thus quanti- 
tatively not as precise as for the harmonic chain, but 
nevertheless yield reliable order of magnitute estimates. 

Let us first analyze condition (|23|l . Since it cannot be 
checked for every state \a) we use the stronger condition 



^-^>J8[A2L 



N, 



(50) 



instead, which implies that 1)231) holds for all states \a). 
We require 1)50)1 to be true for all states with energies in 
the range (|2"7j) . It is hardest to satisfy for E^ = E m i n , 
we thus get the condition on n: 



eo 



(51) 



where e mm = E min /n and e = E /(nN G ). 

We now turn to analyze condition 1)25(1. Equation l|49|) 
shows that the A 2 do not contain terms which are pro- 
portional to Ep. One thus has to determine, when the 
A 2 arc approximately constant which is the case if 



0- 



L f*J max L f 1 -* 



MJ ruin 



< \E U } ~ \E U ] 

L H-imax L 



(52) 



where [x] maa; and [x] m i„ denote the maximal and minimal 
value x takes on in all states \a). As a direct consequence, 
we get 



N «1 



(53) 



which means that temperature is intensive. Defining the 
quantity e M = E^/n, we can rewrite 1)52)1 as a condition 
on n, 



n > 



3 \A 2 } - [A 2 1 

_ L A 1 -I max L M-l mil 

[eJ -\e u ] . 

L h-i max L /' . /;/ in 



(54) 
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where the accuracy parameter <5 1 is equal to the ratio 
of the lhs and the rhs of 1|52|) . 

Since equation i|52[l does not take into account the en- 
ergy range (|27(l . its application needs some further dis- 
cussion. 

If the occupation number of one mode of a group is 
changed, say from = to n^.(u) = 1, the corre- 

sponding differ at most by 4B 2 \K 2 — L 2 \ / (n + 1). 
On the other hand, \A 2 1 - \A 2 ] . = B 2 I K 2 - L 2 1 . 

L A*-l max L PJ mm I I 

The state with the maximal A 2 and the state with the 
minimal thus differ in nearly all occupation num- 
bers and therefore their difference in energy is close to 
i E ^max " \ E n\ m in- On tlie °th er hand, states with simi- 
lar energies also have a similar A^. Hence the A^ only 
change quasi continuously with energy and equation I52|) 
ensures that the A^ are approximately constant even on 
only a part of the possible energy range. 

We are now going to discuss three special coupling 
models. 







K = L = ( 




10 lr 






10 o> 


\ a = 100 


t = l\ 








FIG. 5: Log-log-plot of a function of T/B from eq. 

l|51[l for two values of the accuracy parameter a, a = 1 and 
a = 100, left plot for K = L = 0.1 and right plot for K = 
L — 10. a is denned in equation l|27^ . 



curacy parameter" S. The dependence of the results 
on a is shown in figure [SJ a plays a role only where 
Ermn = E '/ (aN G ) + E /N G (cf. eq. (23). Then 
for smaller a, n m i n eventually decays steeper and thus 
reaches n mm = 1 already at lower temperatures. There 
is thus a temperature interval, where n rnm is larger for 
larger a and vice versa. This dependency has the same 
interpretation as for the harmonic chain. 



A. Coupling with constant width A a : 



If one of the couplings vanishes ( J x 
K = L and A 2 = B 2 K 2 



or J v 



0). 



L and A 2 = B 2 K 2 is constant. In this case only 
criterion (|23() has to be satisfied, which then coincides 

with jni). 

Plugging expressions (glj), (gU and (gHl with J x = J 
and J y = into condition (|51J) . one can now calculate 
the minimal number of systems per group. 





10 8 




10 6 


N ^x 


10 4 


\ 


10 2 



10" 



10" 



10" 



10" 



io- 



T/B 



FIG. 4: Log-log-plot of n m in from eq. 1)51 ft for K = L = 0.1 
(dashed line) and for K = L = 10 (solid line) as a function of 
T/B. ol = 10 is defined in equation l !37t . 

Figure 0] shows n mm for weak coupling K = L = 0.1 
and strong coupling = L = 10 with a = 10 as a 
function of T/B. We choose units where Boltzmann's 
constant fcs is one. 

For any set of parameters, there is a finite temperature 
above which n m i n = 1. 

Note that, since A^ = const, condition (|51J) coincides 
with criterion i|23|) (A M = const = [A^l^J, so that 
using H51f) does not involve any approximations. 

As condition l|24|) is automatically satisfied for the 
present model, the results do not depend on the "ac- 







B. Fully anisotropic coupling: J x — —J y 

If both couplings are nonzero, the variances A^ are 
not constant. As an example, we consider here the fully 
anisotropic coupling, where J x = — J y , i. e. K = 0. Now 
criteria (|51(l and Ij54(l have to be met. 

For K = 0. one has [A^l = B 2 L 2 , \A 2 ] 

L A* J max L AM mm 

and lej = — \e u ] = B. 

L ^imax L MJmin 

Plugging these results into l(54"|) as well as lf4~4"|l and l|4"5"|) 

into H51(l . the minimal number of systems per group can 
be calculated. 

Figure shows n mm from criterion (|51|l and from cri- 
terion l|54(l separately, for weak coupling L = 0.1 and 
strong coupling L = 10 with a = 10 and 5 = 0.01 as 
a function of T/B. For each coupling strength L, the 
stronger condition, that is the higher curve in figure [5] 
sets the relevant lower bound to the group size n. 







- 10 s 






- 10 6 


\ v \ 

N 

\ ^ 




- 10 4 




•n \ 


- 10^\ 



10" 



10" 



10" 



10 2 



io- 



T/B 



FIG. 6: Log-log-plot of n m i n for L = 0.1 from eq. 15111 
(dashed line) and from eq. 1541 (dash - dotted line) and n m ; n 
for L = 10 from eq. 15111 (solid line) and from eq. 1541 (gray 
line) as a function of T/B. K = 0, a = 10 and 8 = 0.01. a 
and 5 axe defined in equations 12711 and 15411 respectively. 
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In the present case, all occupation numbers are 
zero in the ground state of a group. In this state, is 
maximal (A^ = B 2 L 2 ) as can be seen from l|49ll . There- 
fore criterion (|51(l is equivalent to criterion (|23|l for low 
temperatures, where E m i„ = [Efj] min . For high tempera- 
tures, where E m i n = E/(aNc), condition l|51[) is slightly 
stronger than l|23f) . For the present model, this is only 
the case for L = 0.1 (dashed line) and T > 0A5B. 

In figureEl the results obtained from equation (|54*|) arc 
proportional to 8" 1 (dash - dotted line and gray line), 
while those obtained from equation Q51|l (dashed line and 
solid line) have the same dependency on a as shown in 
figure [SJ 

C. Isotropic coupling: J x — J y 





-10 12 




-10 10 




- 10 8 




- 10 6 




MO 4 


N. \ 

^ \ 

^* N. \ 




^ \ 

\ 





1Q -6 1Q -4 1Q -2 1Q 2 1Q 4 



FIG. 7: Log-log-plot of n m i n for K = 0.1 from eq. I|55|l 
(dashed line), and from eq. 1541 (dash - dotted line) and 
fimin for K = 10 from eq. I|51ft (solid line) and from eq. I|54[l 
(gray line) as a function of T/B. L — 0, a = 10 and 8 = 0.01. 
a and 8 are defined in equations 1271 and 1541 respectively. 



As a third example, we consider the isotropic coupling, 
where J x = J y , i. e. L = 0. Again, both criteria H51|) 
and H54f> have to be met. 

The values of \A 2 ] , [A?.] . , [ej and [ej . 

L Wm ni 1 L MJ min > L W mm L Pimm 

are given in equations (|B11|) . (|B12fl and i|B13|l 

For the present model with L = and \K\ < 1 all oc- 
cupation numbers n^(fj.) are zero in the ground state and 
thus A^ = 0. As a consequence, condition l|51|l cannot 
be used instead of lt2~5|) . We therefore argue as follows: 
In the ground state E^ — Eq /Nq — as well as A^ = 
and all occupation numbers are zero. If one occu- 

pation number is then changed from to 1, A^ changes 
at most by 4 B 2 K 2 / (n + 1) and E^ changes at least by 
2 B (l-\K\). Therefore will hold for all states except 
the ground state if 

n>2B(3 l^\K\ (55) 

If l^l > 1, occupation numbers of modes with cos(/c) < 
are zero in the ground state and occupation num- 
bers of modes with cos(fc) > are one. A^ for the 
ground state then is [A?,l ~ lA^l 12 and (f3TI> is a 

° L rJjs L max ' 1 ' 

good approximation of condition (|23|l . 

Plugging these results into (|54|) as well as l|44|) and 
ffSJl into (|5TJ for \K\ > 1 and using for \K\ < 1, the 
minimal number of systems per group can be calculated. 

Figure [7] shows n min from criteria lj5"5|) and (|54"fl for 
weak coupling K = 0.1 and from criteria l|51|) and H54|) 
for strong coupling K = 10 with a = 10 and 8 = 0.01 
as a function of T/B. For each coupling strength K, the 
stronger condition, that is the higher curve in figure \7\ 
sets the relevant lower bound to the group size n. 

Equation (|55(l does not take into account the relevant 
energy range (|27|l . it is therefore possible that a weaker 
condition could be sufficient in that case. However, since 
(|54|l is a stronger condition than I55|) for K = 0.1, this 
possibility has no relevance. 

For strong coupling, K = 10 (|51f) is used to approxi- 
mate (12.'it . This approximation is expected to be good 
because A M is close to its maximal value for low energy 



states. Furthermore, the temperature dependence we ob- 
tain here for n m i n for low temperatures is the same as 
for the harmonic chain, n rn i n oc T~ 3 . This agreement 
is to be expected: The two couplings, when expressed in 
creation and annihilation operators, have the same struc- 
ture and the upper limit of the spectrum of the spin chain 
becomes irrelevant at low temperatures. 

For the present model, the dependence of the results 
on the "accuracy parameters" a and <5 is as follows. Re- 
sults obtained from equation l)54|l are proportional to 
(dash - dotted line and gray line), while the result ob- 
tained from equation (|51() (solid line) has the same de- 
pendency on a as shown in figure [S] For weak coupling 
and low temperatures (dashed line) n min does not de- 
pend on the two "accuracy parameters" . 

V. SUMMARY AND CONCLUSIONS 

We have considered a linear chain of particles interact- 
ing with their nearest neighbors. We have partitioned the 
chain into identical groups of n adjoining particles each. 
Taking the number of such groups to be very large and 
assuming the total system to be in a thermal state with 
temperature T we have found conditions (equations l|23|) 
and H25|) ) , which ensure that each group is approximately 
in a thermal state. Furthermore, we have determined 
when the isolated groups have the same temperature T, 
that is when temperature is intensive. 

The result shows that, in the quantum regime, these 
conditions depend on the temperature T, contrary to the 
classical case. The characteristics of the temperature de- 
pendence are determined by the width A a of the distri- 
bution of the total energy eigenvalues in a product state 
and its dependence on the group energies E a . The low 
temperature behavior, in particular, is related to the fact 
that A Q has a nonzero minimal value. This fact does not 
only appear in the harmonic chain or spin chains but is 
a general feature of quantum systems composed of inter- 
acting particles or subsystems. The commutator [H, Hq] 
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is nonzero and the ground state of the total system is en- 
ergetically lower than the lowest product state, therefore 
A a is nonzero, even at zero temperature [251 IsH Is^. lifij . 

We have then applied the general method to a har- 
monic chain and several types of Ising spin chains. For 
concrete models, the conditions (|2*5|) and determine 
a minimal group size and thus a minimal length scale 
on which temperature may be defined according to the 
temperature concept we adopt. Grains of size below this 
length scale are no more in a thermal state. Thus tem- 
perature measurements with a higher resolution should 
no longer be interpreted in a standard way. 

We have given order of magnitude estimates for the 
minimal group size (minimal length scale) for the models 
mentioned above. The most striking difference between 
the spin chains and the harmonic chain is that the energy 
spectrum of the spin chains is limited, while it is infinite 
for the harmonic chain. 

For spins at very high global temperatures, the to- 
tal density matrix is then almost completely mixed, i. 
e. proportional to the identity matrix, and thus does 
not change under basis transformations. There are thus 
global temperatures which are high enough, so that local 
temperatures exist even for single spins. 

For the harmonic chain, this feature does not appear, 
since the size of the relevant energy range increases indef- 
initely with growing global temperature, leading to the 
constant minimal length scale in the high energy range. 

For the spin chain with isotropic coupling, J x = J y , 
and the harmonic chain, the temperature dependencies 
of Ti m i n for low temperatures coincide, n TO j„ cx T~ 3 , be- 
cause both couplings have the same structure and the 
upper limit of the spectrum of the spin chain becomes 
irrelevant at low temperatures. 

The spin chain with J x = or J y = shows the inter- 
esting feature that A 2 a is constant and condition (|25|1 is 
automatically fulfilled. 

The set of models we have discussed is by no means 
exhaustive. It would be particularly interesting to see 
whether there are systems for which local temperatures 
can exist although they are not intensive. This can hap- 
pen if cither e a or were proportional to E a . A 2 how- 
ever has dimension energy squared, so that it cannot be 
proportional to E a unless there exists another character- 
istic energy of the system independent of E a . So far, we 
have not found models where e a oc E a . 

For the models we consider here, the off diagonal ele- 
ments of the density operator in the product basis, (a\p\b) 
(a 7^ b), are significantly smaller than the diagonal ones, 
(a\p\a). Our general result, conditions (|23[) and (|25[l . 
thus states that the density matrix p "approximately" 
factorizes with respect to the considered partition. This 
implies that the state p is not entangled with respect 
to this partition, at least within the chosen accuracy. 
It would therefore be interesting to see how our result 
relates to the scaling of entanglement in many particle 
systems |38|. 

Unfortunately, our approach only applies to nonzero 



tem per atures. The underlying central limit theorem 
|29l l31| is about the weak convergence of the distribu- 
tion of energy eigenvalues. Weak convergence means that 
only integrals over energy intervals of nonzero length do 
converge. We thus cannot make statements about a sys- 
tem in its groundstate let alone about the entanglement 
in that state. 

Since harmonic lattice models in Dcbye approximation 
have proven to be successful in modeling thermal prop- 
erties of insulators (e.g. heat capacity) j34|, our calcula- 
tion for the harmonic chain provides a first estimate of 
the minimal length scale on which intensive temperatures 
exist in insulating solids, 

Let us give some numerical estimates: Choosing the "ac- 
curacy parameters" to be a = 10 and 6 = 0.01, we get 
for hot iron (T»6k 470 K, a Q m 2.5 A) l mm m 50 pm, 
while for carbon (9 w 2230 K, a w 1.5 A) at room tem- 
perature (270 K) l m i n ~ 10 pm. The coarse-graining will 
experimentally be most relevant at very low tempera- 
tures, where l m in may even become macroscopic. A per- 
tinent example is silicon (0 ks 645 K, ao ~ 2.4 A), which 
has lmin ~ 10cm at T ~ IK (again with a = 10 and 
8 = 0.01). 

Of course the validity of the harmonic lattice model 
will eventually break down at finit, high temperatures 
and our estimates will thus no longer apply there. 

Measurable consequences of the local breakdown of the 
concept of temperature and their implications for future 
nanotcchnology are interesting questions which arise in 
the context of the present discussion. 

In the secnarios of global equilibrium, which we con- 
sider here, a temperature measurement with a micro- 
scopic thermometer, locally in thermal contact with the 
large chain, would not reveal the non existence of local 
temperature. One can model such a measurement with a 
small system, representing the thermometer, coupled to 
a heat bath, representing the chain. It is a known result 
of such system bath models ^lj , that the system always 
relaxes to a thermal state with the global temperature of 
the bath, no matter how local the coupling might be. 

This, however, does not mean that the existence or non 
existence of local temperatures had no physical relevance: 
There are indeed physical properties, which are deter- 
mined by the local states rather than the global ones. 
Whether these properties arc of thermal character de- 
pends on the existence of local temperatures. A detailed 
discussion of such properties will be given elsewhere. 

The length scales, calculated in this paper, should also 
constrain the way one can meaningfully define tempera- 
ture profiles in non-equilibrium scenarios |42| . Here, tem- 
perature measurements with a microscopic thermome- 
ter, which is locally in thermal contact with the sample, 
might indeed be suitable to measure the local temper- 
ature. An explicit study of this possibility should be 
subject of future research. 
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APPENDIX A: DIAGONALIZATION OF THE 
HARMONIC CHAIN 

The Hamiltonian of a harmonic chain is diagonalized 
by a Fourier transformation and the definition of creation 
and annihilation operators. 

For the entire chain with periodic boundary conditions, 
the Fourier transformation reads 



1j 
Vj 



1 «t exp(ia kj) 

y/nNc k \ v k exp(-ia kj) 



(Al) 



with k = 2Trl/(a nN G ) and (I = 0,±1,..., 
±(niV(3 — 2)/2, (nN G )/2, where nN G has been assumed 
to be even. 

For the diagonalization of one single group, the Fourier 
transformation is 



1i I = 
Pj J 



sin(a fcj) 



(A2) 



with k = nl I (ao (n + 1)) and (Z = 1,2,..., n). 

The definition of the creation and annihilation opera- 
tors is in both cases 



"A- 



k K = 



1 



\/2rruJ k 



muj k u k 



iv k 



(A3) 



where the corresponding u k and v k have to be inserted. 
The frequencies uj k are given by Lo k = 4wg sin 2 (fcao/2) in 
both cases.. 

The operators a k and a k satisfy bosonic commutation 
relations 

[ak,a p ] = 
[a k ,al] = S kp (A4) 

and the diagonalized Hamiltonian reads 



H = J2oj k (a\,a k + ]- 



(A5) 



APPENDIX B: DIAGONALIZATION OF THE 
ISING CHAIN 

The Hamiltonian of the Ising chain is diagonalized 
via Jordan- Wigner transformation which maps it to a 
fermionic system |3^, IU 




(Bl) 



The operators c; and c, 
commutation relations 



fulfill fermionic anti- 



{d,cj} = 
{c u c\} = 6 i: 



(B2) 



and the Hamiltonian reads 



H 



D 



E 



2c]c, 



3 



c]c j+1 + h.c. 



c+ct 

C 3 C 3+1 



h.c 



(B3) 



with K = (J x + J y )/(2B) and L = (J x — J y )/(2B). 
In the case of periodic boundary conditions a bound- 
ary term is neglected in equation (|B3|) . For long chains 
(nNc — ► 00) this term is suppressed by a factor (nNc)~ 1 . 
The Hamiltonian now describes Fermions which interact 
with their nearest neighbors. As for the bosonic system, a 
Fourier transformations maps the system to noninteract- 
ing fermions. For the whole chain with periodic boundary 
conditions 



1 



G 



4 



(B4) 



with k = (2nl)/(nN G ) where I = 0,±1, 
±(nN G - 2)/2, {nN G )/2 for nN G even, and 



- ^ sin(fcj) 



4 



(B5) 



with k = (nl)/(n + 1) and (I = 1, 2, . . . , n) for one single 
group. 

In the case of periodic boundary conditions, fermion 
interactions of the form d\d}_ k and d k d^ k remain. There- 
fore, one still has to apply a Bogoliubov transformation 
to diagonalize the system, i.e. 



u k b\ - Wkb-k 



d k = u k b k + iv k W_ k 



where u k 



U-k, v k 



-v- k 



and 



(B6) 



1. With 



the definitions u k = cos(Q k /2) and v k = sin(9/c/2) the 
interaction terms disappear for 



cos(e fe ) 



1 — K cos k 



- A' cos k} 2 + [L sinfc]' 



(B7) 



In the case of the finite chain of one group, the Bogoli- 
ubov transformation is not needed since the correspond- 
ing terms are of the form d k d k and d k d k and vanish by 
virtue of equation (|B2|) . 

The Hamiltonians in the diagonal form read 



1 



(B8) 
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where the frequencies are 

uj k = 2B v/[l - K cos k] 2 + [L sin k} 2 (B9) 

with k = (2itl) / (tiNg) for the periodic chain and 

uj k = 2B(1 - Kcosk) (BIO) 

with k = (irl)/(n + 1) for the finite chain. 

For the finite chain the occupation number operators 
may also be chosen such that u) k is always positive. Here, 
the convention at hand is more convenient, since the same 
occupation numbers also appear in the group interaction 
and thus in A^. 

1. Maxima and minima of _E M and 

The maximal and minimal values of are given by 

fA max ) = \t)nB, (Bll) 



for I If I < 1 and by 



\E, 



nB- 



Vi< 2 - 1 



arcsm 



(B12) 



for | If | > 1, where the sum over all modes k has been 
approximated with an integral. 

The maximal and minimal values of A 2 are given by 



l H'-imax 

A 2 

L w, 



= B 2 



max (K 2 , L 2 ) 
min (K 2 ,L 2 ) 



(B13) 
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